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The definition of a cluster state naturally suggests an implementation scheme: find a physical 
system with an Ising coupling topology identical to that of the target state, and evolve freely for a 
time of Yj. Using the tools of optimal control theory, we address the question of whether or not 
this implementation is time-optimal. We present some examples where it is not and provide an 
explanation in terms of geodesies on the Bloch-sphere. 
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I. INTRODUCTION 



The one-way quantum computation model is an ap- 
proach to quantum information processing where the evo- 
lution is driven by local operations and measurements 
only P, In the experimental realization of such a 
model, preparation of the highly entangled initial state 
is therefore of primary concern. These initial states - so 
called cluster states - are also interesting in their own 
right, due to the favorable scaling of their entanglement 
properties [1, To date, cluster states of four to six 
qubits have been realized experimentally in photonic sys- 
tems 0, and are also actively pursued in other ar- 
chitechtures, e.g. ion traps 

In this paper, the efficient generation of cluster states 
is studied using techniques from optimal control theory 
Given an experimental framework, the aim of these 
techniques is to find the optimal set of controls to steer 
the system so that a desired target state or unitary gate 
is implemented. For systems consisting of two qubits, 
general analytical solutions exist for the construction of 
time-optimal unitary transformations [l^, [HI and state- 
to-state transfers [12, [ll] if fast local controls are avail- 
able. For three or more qubits, analytical solutions are 
only known in some special cases, it was 

shown that the time-optimal generation of indirect cou- 
plings and trilinear Hamiltonians and the efficient trans- 
fer of order along Ising spin chains [H, [l^, [l^l can be 
reduced to the problem of computing singular geodesies. 
In addition, powerful numerical methods [2l| are avail- 
able that make it possible to explore the physical limits 
of time-optimal control experiments in cases where no 
analytical solutions are known. The GRAPE algorithm 



21 1 has been used for the generation of quantum gates 
22 1 and state-to-state transfers in NMR [23|, as well as 
for superconducting qubits [22]. Here, analytical and nu- 
merical techniques are used to investigate the problem of 
time-optimal cluster state preparation. 



A cluster state is defined by a graph. To prepare the 
n-qubit cluster state corresponding to a graph G: 

1. Prepare (locally) the initial state 



1^) := 



+ 11) 



V2 



2. Evolve under the Ising Hamiltonian 
7rJ(i) 



H. = 



Y: + (1) 



(a, a') 



for a time such that J{t)dt = i. The sum is 
over all edges of G, where each edge connects the 
qubit pair (a, a'). 

We consider 7i-qubit systems of the form 

J 

where the Uj are time-dependent functions to be chosen, 

and the He characterize the available controls. Two 
different control settings will be considered: 

(i) Local X and y control on each qubit: 

Hi'^''^ = la(j\ = (2) 

where j G {1, 2, 2n}. 

(ii) A single global x control: 



U) 



(3) 
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In the following, we allow for fast local controls, i.e. the 
functions Uj are unrestricted. 
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II. 3 COMPLETELY-COUPLED QUBITS 

To shed some light on when a speedup may be pos- 
sible, we begin with a symmetric 3-qubit system which 
will prove analytically tractable. The qubits are Ising- 
coupled according to the complete coupling graph K3 
[2^ . illustrated in fig. [T] Without loss of generality we 



approximately 0.77 x -j, as seen in fig. [2^. The signif- 
icance of this particular value will become clearer later. 
Fig. [2)3 shows a numerically optimised solution close to 
the minimal time. We find in all observed cases that 
the maximum fidelity does not depend on whether the 
controls are specified according to ([2|) or 

Motivated by these symmetric solutions, we now re- 
strict ourselves to control setting (ii), specifying He ac- 
cording to ([3]). This control Hamiltonian, in addition to 
the drift Hamiltonian and the inital and target states, is 
symmetric under the cyclic permutation operator 



FIG. 1: coupling graph. 

will drop the local terms in ([1]), as only the entangling 
part of the operation contributes to the time required. 
Furthermore J is assumed to be constant. The clus- 
ter state is then defined as 

in) ■■= exp (-^^H^ \h) , 



where 



ttJ 



The preparation of this state poses the following control 
problem: maximize 



F{U) :=|(r3|.(C/|/3))| 



(4) 



subject to equation of motion U = —iHU. In the first 
instance we will allow for full local control on the qubits 
and thus specify He according to ([2]). 

To maximize the fidelity defined in ^ in the short- 
est possible time, we first appl y a numerical gradient- 
ascent algorithm, as detailed in [2l| . Over a comprehen- 
sive range of initial conditions the controls are updated 
iteratively, incrementing _F to a local maximum. The 
minimal time found by the algorithm to generate IT3) is 




2Jt 
(b) 



FIG. 2: (a) Maximum achievable fidelity as a function of the 
transfer time, (b) A sample solution at the minimal time 
for initial controls = Uy{t) = 0. The y controls are 

omitted as they remain zero, and the x controls are identical 
on each qubit due to the permutation symmetry of the drift 
Hamiltonian and the initial and target states. 
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and the persymmetry operator 

P ■■= i<J.f' , 

which themselves commute. The dynamics are thus re- 
stricted to the simultaneous eigenspace of S and P cor- 
responding to the eigenvalue pair {1,1}. Transforming 
to a new basis composed of the simultaneous eigenstates 
of S and P makes this explicit [2^. The Hamiltonians 
are now diagonalized into 2x2 blocks. As we need only 
consider the {1,1} block, the state transfer problem can 
be reduced to 



14) 



V3 
1 



V3 

-1 



under the Hamiltonian H' = H'^^ + u{t)H'e, where 



H',^^-l 



-1 
3 



" 2 



2 V3 
\/3 



(5) 



(6) 



In order to consider the problem geometrically, we rep- 
resent the transfer on the Bloch sphere by projecting onto 
the axes Ij :— (Jj/2, as illustrated in fig. [3l The initial 
and final states j/g) and \T^ are identified with vectors 
a = ^(\/3, 0, 1) and b = \/3, 0, 1), respectively, while 
the Hamiltonians i?^ and H'^ correspond to (unnormal- 
ized) rotation axes (0,0, — 27rJ) and (-\/3, 0,1), respec- 
tively. The trivial solution is to set u{t) — Q and rotate 
about H'j^ for ^ units of time, but faster solutions may 
exist. Motivated by the numerical results in Fig. ^Bp, we 
first restrict u{t) to the following form: 

1. Constant pulse u over time interval [0,T]. 

2. Hard pulse of angle 4> at time T. 

After deriving a time-optimal solution in this setting, we 
will show that it remains time-optimal when the restric- 
tions are removed and general time-varying pulses are 
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(0,0,1) 



C=(-V3/2,0,1/2) 




terms of 9 as 



a = (V3/2,0,1/2) 



FIG. 3: (Color online) Transfer of eqn. ([5]) on the Bloch 
sphere. The u = solution (blue line) transfers a to c in a 
time of A time-optimal solution (red line) transfers a to 
h in , followed by a hard pulse (dashed red line) from h 
to c in negligible time. 



considered. The solution is 



(7) 



which explains the value of T « 0.77 x jj obtained nu- 
merically. To demonstrate the time-optimality of this 
solution, we can consider the geometric constructions in 
fig. m Starting from a, the task is to transfer the state to 
any point h on the circle bee obtained by rotating about 
H'^. The hard pulse then transfers 5 to c in an arbitrarily 
small time. The choice of constant u specifies a rotation 
axis H\ at an angle 9 to H'^. The transfer time to be 
minimized is 

„ Zad6 



Zad& = 2 arcsin 



so that the time is 



2 sin 6* 



sm( 



sm a arcsin 



V3 ^ 
2 sin 6* 



Noting that 9 must lie in the interval [-1, ^] for inter- 
section with the circle, we find the maximum at 9 = ^. 
The time required is Tmin — ■ "^^'^ angle of the final 
hard pulse is — -I in this picture, but since \H'^\ — 2 the 
angle around a normed axis is (/> = — f . 

It remains to show that bee cannot be reached in less 
than Tinin when allowing for time- varying pulses. For this 
we introduce 



\H' 



which is simply the rotation axis orthogonal to H'^. We 
consider a generic time- varying control u{t). Our aim 
is to compare each segment of this generic path to the 
optimal one. Let ci and C2 be two circles generated by 
rotating about H'^, chosen to be close enough so that 
u{t) is well approximated by a constant in the interval 
between them. We consider the time required to travel 
from ci to C2 along two different paths: our proposed 
optimal solution ai 02, obtained by rotating purely 
about -ff"*", and a generic path 61 — > 62- Suppose the 
evolution along 61 — > 62 takes a time A. This evolution 
can be decomposed according to the Trotter formula 



where Aadb is the angle swept out by the Bloch vector 
and the length \H'\ gives its angular velocity of rotation. 
Using simple geometry these quantities are expressed in 



(8) 



which is represented graphically in fig. [5^. Note that the 
time required for the operation on the right-hand side of 




FIG. 4: Geometric constructions used. All points shown here 
lie in the Ix-Iz plane except for b, which lies above it on the 
upper surface of the sphere. 




FIG. 5: (Color online) (a) An arbitrary path can be decom- 
posed into rotation about H'^ and rotation about the orthogo- 
nal axis via the Trotter decomposition, (b) Rotating our 
viewpoint by 90°, we see that the optimal trajectory from 
tsi a2 minimizes the angle rotated through in each seg- 
ment (red), when compared to a generic trajectory (blue), ie. 

01 < 6*2. 
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([8]) is still A, as the evolutions along H'^ can be arbitrarily 
fast. The time in a single segment is equal to the angle 
swept out divided by the norm of H^. This angle is 
minimized in every segment when travelling from ai — > 
02, as fig. [SJd illustrates. The time-optimal solution is 
therefore to rotate purely about H^, which corresponds 
exactly to solution ^ . The speedup here arises from the 
fact that H'^ and H'^ are not orthogonal, which can be 
carried over into higher dimensions. 

III. HIGHER DIMENSIONAL CASES 

Finally we provide the minimal times for cluster state 
preparation on a variety of graphs, some of which are il- 
lustrated in fig. [S] to clarify our notation. The times are 



• • • • • • 




• • • • • • 

(a) (b) (c) 



sists between the non-orthogonality of H'^ and H'^ and 
the minimal time - in all cases considered, a speedup is 
possible when H'^ and H'^ are non-orthogonal. Minimal 
time solutions found by the algorithm for Kc, and ify are 
shown in fig. [71 We observe that the effective hard rota- 




(a) (b) 

FIG. 7: Numerically optimised sequences for (a) Ks and (b) 
Kr. 

tion pulse at time T is still present, while the preceding 
pulse shape is no longer constant. 



FIG. 6: Some of the graphs considered: (a) C4, (b) K4, (c) 

G2,3- 

evaluated numerically using the GRAPE algorithm, and 
are included in Table |T1 For all of the 4-qubit graphs, 
both control schemes (i) and (ii) were considered, yield- 
ing exactly the same minimal times in each case. For 
the larger graphs only control scheme (ii) was consid- 
ered, allowing us to reduce the problem to a dimension d 
using a symmetry- adapted basis, as in Sec. HD The times 
listed here hold not just for the target state but its entire 
local unitary orbit, which may include other entangled 
states of interest [23|. While a Bloch sphere analysis is 
not possible for d > 2, we find that a connection per- 



Graph d \{H'^\H'c)\ Tmin 
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0.2 


0.77 


L3 


3 





1.00 


K4 


3 





0.91 


C4 


4 





1.00 


Kb 


3 


0.1 


0.70 


Ke 


4 





1.00 


G2,3 


14 





1.00 


Kr 


10 


0.07 


0.60 



IV. SUMMARY 

For several known two-qubit quantum gates, the 'do 
nothing' operation, i.e. the evolution under a given cou- 
pling Hamiltonian, is time-optimal. For example, this is 
the case for a SWAP gate in the presence of an isotropic 
Heisenberg coupling [28|. Hence it may be surprising to 
find that it is possible to create cluster states faster than 
the time required by the straightforward implementation 
their definition implies. Here we provide examples where 
this is the case. In particular for three Ising-coupled 
qubits with identical coupling constants, the problem of 
time-optimal cluster state generation could be solved an- 
alytically via geodesies on a sphere. The numerical tech- 
niques used here can also be applied to the problem of 
cluster state generation in an arbitrary coupling topology, 
where real experimental values for the coupling constants 
can be used. This will be addressed in a future work. 



Acknowledgements 



TABLE I: Minimal times calculated by the GRAPE algo- 
rithm, with an estimated numerical accuracy of ±1 on the 
last digit. As for the 3-qubit example under control setting 
(ii), included in the first row as a reference, the drift and con- 
trol can be reduced to matrices H'^ and H'^ of size d x d. H' 
indicates a rescaling to unit norm. 



This work has been supported by the Elite Network 
of Bavaria programme QCCC, the integrated EU project 
QAP, and the Deutsche Forschungsgemeinschaft incen- 
tive SFB-63I. R.F. and S.G. would like to thank Christof 
Wunderlich and Thomas Schulte-Herbriiggen for valuable 
discussions. 



5 



[1] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 
5188 (2001). 

[2] M. A. Nielsen, e-print arXi v:quant-ph/0504097[ Rev. 
Math. Phys. (to be pubhshed). 

[3] H. J. Briegel and R. Raussendorf, Phys. Rev. Lett. 86, 
910 (2001). 

[4] W. Diir and H. J. Briegel, Phys. Rev. Lett. 92, 180403 

(2004) . 

[5] P. Walther, K. J. Resch, T. Rudolph, E. Schenck, H. We- 

infurter, V. Vedral, M. Aspelmeyer, A. Zeilinger, Nature 

434, 169 (2005). 
[6] C.-Y. Lu, X.-Q. Zhou, O. Giihne, W.-B. Gao, J. Zhang, 

Z.-S. Yuan, A. Goebel, T. Yang, J.-W. Pan, Nature 

Physics 3, 91 (2007). 
[7] H. Wunderlich and C. Wunderlich, e-print 

[arXiv:quant-ph/0901088 1 . 
[8] L. Pontryagin, B. Boltyanskii, R. Gamkrelidze, E. Mish- 

chenko. The Mathematical Theory of Optimal Processes, 

Wiley-Interscience, New York (1962). 
[9] A. E. Bryson and Y.-C. Ho, Applied Optimal Control, 

Hemisphere, Washington, D.C. (1975). 
[10] N. Khaneja, R. Brockett, S. J. Glaser, Phys. Rev. A 63, 

032308 (2001). 

[11] H. Yuan and N. Khaneja, Phys. Rev. A 72, 040301(R) 

(2005) . 

[12] T. O. Reiss, N. Khaneja, S. J. Glaser, J. Magn. Reson. 

154, 192 (2002). 
[13] N. Khaneja, F. Kramer, S. J. Glaser, J. Magn. Reson. 

173, 116 (2005). 
[14] N. Khaneja, S. J. Glaser, R. Brockett, Phys. Rev. A 65, 

032301(2002). 



[15] T. O. Reiss, N. Khaneja, S. J. Glaser, J. Magn. Reson. 
165, 95 (2003). 

[16] N. Khaneja, B. Heitmann, A. Sporl, H. Yuan, T. Schulte- 
Herbriiggen, S. J. Glaser, Phys. Rev. A 75, 012322 

(2007) . 

[17] H. Yuan, S. J. Glaser, N. Khaneja, Phys. Rev. A 76, 
012316 (2007). 

[18] H. Yuan, R. Zeier, N. Khaneja, Phys. Rev. A 77, 032340 

(2008) . 

[19] N. Khaneja, S. J. Glaser, Phys. Rev. A 66, 060301(R) 
(2002). 

[20] D. Stefanatos, S. J. Glaser, N. Khaneja, Phys. Rev. A 

72, 062320 (2005). 
[21] N. Khaneja, T. O. Reiss, C. Kehlet, T. Schulte-Herb- 

riiggen, S. J. Glaser, J. Magn. Reson. 172, 296 (2005). 
[22] T. Schulte-Herbriiggen, A. Sporl, N. Khaneja, S. J. 

Glaser, Phys. Rev. A 72, 042331 (2005). 
[23] J. L. Neves, B. Heitmann, T. O. Reiss, H. H. R. Schor, 

N. Khaneja, S. J. Glaser, J. Magn. Reson. 181, 126 

(2006). 

[24] A. Sporl, T. Schulte-Herbriiggen, S. J. Glaser, V. Berg- 
holm, M. J. Storcz, J. Ferber, F. K. Wilhelm, Phys. Rev. 
A 75, 012302 (2007). 

[25] D. B. West, Introduction to Graph Theory, Prentice Hall, 
Upper Saddle River, NJ, 2001. 

[26] C. N. Banwell and H. Primas, Mol. Phys. 6 C, 225 (1963). 

[27] M. Hein, J. Eisert, H. J. Briegel, Phys. Rev. A 69, 062311 
(2004). 

[28] T. Schulte-Herbriiggen and O. S0rensen, Concepts Magn. 
Reson. 12, 389 (2000). 



